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ABSTRACT 


The  spectral  profile  for  the  absorption  of  infrared  laser  light  by  a 
vibrational  bond  between  a  physisorbed  atom  and  a  harmonic  crystal  is  calculated. 
We  obtained  an  analytical  expression  for  the  line  shape,  which  includes  the  finite 
memory- time  effects  in  the  interaction  between  atomic  motion  and  bulk- atom 
vibrations.  Both  the  memory  in  the  time  regression  of  the  dipole  correlation 
function  and  the  initial  correlations  are  taken  into  account.  It  is  shown  that 
absorption  from  a  laser  with  a  frequency  which  is  larger  than  the  cutoff  frequency 
Wp  of  the  dispersion  relation  of  the  crystal  can  only  occur  due  to  a  memory  in  the 
relaxation  process,  provided  that  multiphonon  transitions  are  negligible.  We 
predict  a  resonance -like  line  at  (with  the  unperturbed  resonance)  for 

atom- surface  bonds  with  a  permanent  dipole  moment. 


If  a  crystal  is  exposed  to  an  atomic  vapor,  then  many  atoms  will  stick  to  its 


surface  due  to  the  van  der  Waals  interaction.  Every  atom  in  the  vicinity  of  the 

solid  experiences  an  attractive  potential,  which  supports  continuum  (desorbing) 

states  and  bound  (adsorbing)  states.  Transitions  from  bound  to  continuum  states 

can  be  induced  by  thermal  coupling  of  the  vibrational  bond  with  the  phonon 

reservoir  of  the  crystal,  or  by  illumination  with  a  strong  infrared  laser.  With 

the  first  mechanism  a  large  desorption  rate  can  be  achieved  by  a  sufficient 

heating  of  the  substrate,  whereas  in  the  second  process  the  transition  to  the 

continuum  is  brought  about  by  photon  absorption.  For  these  processes  the 

interaction  between  adsorbate  and  crystal  can  be  accounted  for  by  relaxation  terms 

in  a  master  equation  for  the  populations  of  the  vibrational  levels,  where  the  rate 

constants  are  given  by  the  Golden  Rule.  Dynamical  properties  of  the  system  are 

determined  by  the  relative  values  of  these  rate  constants  and  by  the  level 

structure  of  the  potential.  Then  the  desorption  rate  as  a  function  of  time 

contains  information  on  the  atom-crystal  interaction,  e.g.  the  rate  constants. 

1-3 

Several  authors  applied  this  technique  to  evaluate  the  photodesorption  yield  as 
a  function  of  the  laser  frequency.  A  far  more  sensitive  method  to  obtain  insight 

i 

into  the  details  of  the  dynamical  features  of  adsorbates  is  by  measuring  the 
steady-state  low- intensity  absorption  profile  I(w)  as  a  function  of  the  probe 

(laser)  frequency  w.**  ^  Then  the  absorption  spectrum  will  reveal  the  details  of 

7  8 

dynamical  atom-lattice  bonds,  surface-modified  internal  molecular  modes,  or 

9 

properties  of  the  interaction  between  two  adspecies.  We  shall  consider  a  single 
atom  which  is  adsorbed  on  a  harmonic  -  lattice  crystal,  and  bounded  to  the  surface 
by  a  potential  V(z).  We  neglect  lateral  motion  and  indicate  the  normal  to  the 
surface  by  the  z-axis.  Then  the  induced  dipole  moment  jj(t),  where  t  denotes  the 
Heisenberg  picture,  of  the  atom-crystal  bond  must  be  in  the  z -direction,  and  hence 


we  can  write  ji(t)  -  p(t)ez  with  p(t)  a  scalar  operator.  An  infrared  laser  with 
intensity  IT  (energy  per  unit  of  time  which  passes  a  unit  surface  area, 
perpendicular  to  the  direction  of  propagation)  and  polarization  £^  is  incident  on 
the  atomic  bond.  Since  the  wavelength  of  the  radiation  is  much  larger  than  the 
atom-surface  separation,  we  can  adopt  the  dipole  approximation  for  the  interaction 


between  the  bond  and  the  external  field.  Then  a  general  expression  for  the 
absorbed  energy  per  unit  of  time,  the  absorption  profile,  reads 


I(u)  -  IL|ez-£L| 


— “•  Re  lim  f° 

£o*C  t-~Jo 


dr  elti,rTrp ( t)  [n(t+r)  ,n(t)  } 


(l.D 


where  p(t)  is  the  density  operator  for  the  atom,  the  entire  crystal,  and  the 
interaction.  The  appearance  of  the  commutator  reflects  that  the  net  absorption  is 
a  balance  between  stimulated  absorption  (the  term  /x(t+r  )/*(  t) )  and  stimulated 
emission  (the  term  /i(t)/*(t+r ) )  of  photons  from  and  into  the  laser  field. 

Transformation  of  Eq.  (1.1)  to  the  Schrodinger  picture  yields  the  equivalent 
expression 


I(m)  -  ILl£z-£L|2  7^  Re  lim  dr  TrMel(w'L)r 


M,p(t) 


(1.2) 


Here,  the  limit  t-«>  only  pertains  to  the  density  operator  p(t),  and  obviously 
p( t*-»)  represents  the  thermal-equilibrium  state  of  the  system.  If  we  write  a  for 
an  arbitrary  density  operator,  as  we  shall  do  throughout  the  paper,  then  the 
Liouvillian  L  in  the  exponential  in  Eq.  (1.2)  is  related  to  the  Hamiltonian  H  by 


Lo  -  H  [H,o] 


(1.3) 


From  Eq.  (1.2)  we  notice  that  the  frequency  dependence  of  I(w)  is  governed  by  the 
Fourier- Laplace  parameter  u  in  the  transformation  of  exp(-iLr),  and  consequently 
the  resonance  lines  in  the  profile  are  situated  at  the  eigenvalues  of  L.  Then  we 
recall  that  the  equation  of  motion  for  the  density  operator  is 


1  at  -  Lf 


(1.4) 


and  therefore  the  dynamical  properties  of  the  system  will  be  displayed  in  the  in¬ 
dependence  of  I(w),  even  though  the  system  is  in  a  stationary  state  (for  which 
p(t)  becomes  independent  of  time).  It  is  not  the  direct  time  evolution  of  p(t) 
which  is  probed  by  the  laser,  but  the  time- regression  operator  exp(-iLr)  of  the 
dipole  correlation  functions,  as  they  appear  in  Eq.  (1.1).  The  significance  of  a 
measurement  of  I(w)  then  relies  on  the  fact  that  L  represents  the  entire  system. 


(2.1) 


in  terms  of  the  annihilation  (a^s)  and  creation  (a£s)  operators  for  phonons  in  the 

mode  ks.  Here,  k  and  s  denote  the  wave  vectors  and  polarizations,  respectively, 
and  «s(k)  is  the  dispersion  relation.  The  Hamiltonian  for  the  bounded  atom 
includes  a  kinetic  energy  and  a  potential 
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H 

a 


-2  £♦»<■> 


(2.2) 


with  m  the  mass  of  the  atom.  Eigenstates  and  eigenvalues  of  can  be  found 
easily  for  a  variety  of  potentials  V(z) .  An  important  example  is  the  Morse 
potential,  which  models  the  atom-crystal  binding  quite  accurately.  Kinetic 
coupling  between  the  atomic  motion  and  the  phonon  field  is  assumed  to  be  dominated 
by  single-phonon  interactions,  for  which  the  Hamiltonian  reads 


H 

ap 


v  dv 

-iu*e  )  , 

-  -z  dz 


(2.3) 


The  operator  u  is  the  displacement  of  the  crystal  atom  which  is  closest  to  the 
adsorbate.  Explicitly^ 


ks 


Mv 


k2MVws(k)' 


(aks+a-ks>2ks 


( 2.1) 


with  V  and  v  the  volumes  of  the  crystal  and  a  unit  cell,  respectively,  M  the  mass 
of  a  crystal  atom,  and  the  unit  polarization  vector  of  a  phonon  in  the  mode 
ks.  The  total  Hamiltonian  then  becomes 


H  -  H  +  H  +  H 
p  a  ap 


(2.5) 


which  determines  the  Liouvillian  L  according  to  Eq.  (1.3).  Then  the  density 
operator  p( t)  of  the  system  follows  after  solution  of  Eq.  (1.4),  and  the  spectral 
profile  is  obtained  in  evaluating  expression  (1.2). 


•mww wrwrw yqr^w ur<tw w  tv»v  wir^v  *  w www^gF"gynpngr«*r» w  l  »  u-wy^ w »  jrrjrw 


III.  RESERVOIR  CORRELATION  PJNCTION 

Due  to  the  many  degrees  of  freedom  of  the  phonon  field  a  direct 
diagonalization  of  the  Hamiltonian  is  intractable,  so  that  we  have  to  resort  to  an 
approximation.  The  thermal -equilibrium  density  operator  of  the  crystal  at 
temperature  T  is 


-1 


p  -  {Tr  exp( -H  /kgT) )  exp(-H  /kgT) 


(3.1) 


with  kn  Boltzmann's  constant,  and  where  Tr  indicates  a  trace  over  the  states  of 
B  p 

the  phonon  field  only.  Then  the  idea  is  that  the  large  crystal  can  be  regarded  as 

a  thermal  reservoir,  and  that  its  state  p  is  not  affected  by  the  presence  of  the 

P  .  12 
single  atom  on  the  surface.  The  central  quantity  in  standard  relaxation  theory 

is  the  reservoir  correlation  function 


.9  iLnr  - 

f(r)  -  Y.  Tr  (u*e  )e  P  (p  u*e  )  , 


(3.2) 


with  L  a  —  [H  ,a ]/W.  Relaxation  constants  are  then  expressed  in  the  Fourier- 
P  1  P 

Laplace  transform  of  f (r) 


?( 


«)  -  r 

Jo 


dr  e^wrf(r) 


(3.3) 


We  shall  adopt  a  Debye  model  for  the  dispersion  relation,  which  implies 


cs<k)  -  c'k  H(wD-c'k)  , 


(3.4) 


in  terms  of  the  Debye  frequency  the  speed  of  sound  c' ,  and  the  unit-step 
function  H.  Furthermore,  we  recall  the  relation 


v  6ir 

,3  “  3  ’ 

C  WD 


(3.5) 


expressing  that  the  cutoff  frequency  appears  as  a  consequence  of  the  finiteness 
of  the  volume  v  of  a  unit  cell.  Then  it  is  an  easy  matter  to  compute  f(ci>),  and  we 
obtain^ 


?(w)  -  S  g(w/a>D) 


(3.6) 


where  the  overall  factor  is  given  by 


(3.7) 


The  dimensionless  function  g(w/w^)  which  represents  the  u>- dependence  of  ?(w)  is 


found  to  be 


g(2)  -  zH(z)H(l-z)  -  i  (l+zlog|l-J|)  +  — f—  H  ( 1  - 1  z  |  ) 

*  Z  e7  Z  -1 


i  P  [  dx  — - —  x 
*  JO  e7X-l 


2  2  ’ 
z  -x. 


(3.8) 


which  depends  parametrically  on  the  temperature  through 


7-kBT  ' 


(3.9) 


N- sv-V-'.V, 


I 


£ 


It  will  appear  that  the  function  g(z)  incorporates  the  refinements  of  an 
absorption  line  in  comparison  with  the  standard  Lorentzian.  A  plot  of  g(z)  can  be 
found  in  Ref.  13. 

Two  important  properties  of  g(z),  which  can  be  deduced  from  Eq.  (3.8),  are 


Re  g(z)  -  0 


for  |zj  >  1  , 


(3.10) 


g(-z)  -  g(z)*  for  T  -» 


(3.11) 


Furthermore,  g(z)  goes  to  zero  very  fast  for  |z[  >  1.  This  implies,  in  view  of 
Eq.  (3.3),  that  f(r)  has  a  time  width  of  the  order  of  l/u>_.  In  a  previous  paper^ 
we  pointed  out  that  this  feature  prohibits  the  application  of  the  Markov 
approximation  in  the  derivation  of  an  equation  for  the  reduced  adsorbate  density 
operator  p  (t),  defined  by 


P  (t)  -  Tr  p(t)  . 

a  p 


(3.12) 


If  only  a  master  equation  for  the  populations  of  the  vibrational  states  is  of 
interest  (as  for  instance  in  the  desorption  problem) ,  then  the  finite  time  width 
of  f(r)  does  not  have  much  significance,  but  for  the  evaluation  of  an  absorption 
profile  it  is  of  paramount  importance  that  the  details  of  ?(w)  are  taken  into 
consideration,  as  we  shall  show  below. 


Finite  memory- time  reservoir  theory  is  a  complicated  mathematical  tool,  which 
can  be  applied  to  solve  Eq.  (1.4)  for  p  (t)  and  to  evaluate  steady-state  quantum 

el 

correlation  functions,  as  they  appear  in  Eq.  (1.1).  Recently  we  developed  the 


iwu»n wi  VI  «■’ <.■  \  t~»  7  »  t  ’r-.'  tv  vt  v—  v-,r  ->  ■ry?’-.. 1  j»  -y  ■*  - .« 


I 

£ 


r<* 


general  theory,^  and  in  this  paper  we  apply  the  formalism  to  the  computation  of  j 
line  shapes  of  adsorbates. 

If  we  would  be  able  to  prepare  the  adsorbate  at  time  zero  in  state  p  (0) , 

ci 

then  its  state  for  t  >  0  is  given  by 


pa(w)  “  w-L  +if(w)  ^a(0) 


(4.1) 


in  the  Fourier-Laplace  domain.  As  usual,  L  denotes  the  commutator  with  H 

a  a 

divided  by  h-  Coupling  to  the  reservoir  is  embodied  in  the  relaxation  operator 
T (w) ,  defined  by 


r(w)o  -  Tr  L  — r1—  .  „  v„  y  , 

a  p  ap  w-L  -L  ap  a  p 
a  p 


L  (a  p  )  , 

an  '  r  t-i  /  » 


(4.2) 


where  is  an  arbitrary  adsorbate  density  operator,  and  L  p  “  [Hap,p]/h-  The 
operator  inversion  on  the  right-hand  side  of  Eq.  (4.2)  might  seem  awkward,  but  in 
Ref.  15  we  have  shown  how  to  evaluate  explicitly  the  matrix  elements  of  T(w) . 

For  the  steady- state  line  profile  we  only  need  the  long-time  solution  of 
Pa(t),  which  can  be  found  from  its  Fourier-Laplace  transform  according  to 


Pa  “  Tint  Pa(t)  -  lim  -iwp  =  (u) 


(4.3) 


t-*“> 


W>0 


With  Eq.  (4.1)  we  then  see  that  p  is  the  solution  of 

a 


(L  -  ir(0))p  -  0  , 

*  d 


and  of  course  the  restictions  p^  -  p  ,  Tr  p  -  1  should  be  imposed. 

3  cl  3.  3. 


(4.4) 


] 


In  the  Markov  or  zero  memory-time  approximation,  the  relaxation  operator  F(u) 
acquires  a  frequency- independent  value.  It  is  the  w-dependence  of  T(w)  which 
reflects  the  memory  in  the  atom-crystal  interaction.  In  the  long-time  limit  the 
density  operator  is  determined  by  T(a>)  at  u  -  0,  but  this  operator  is  not  equal  to 
its  Markovian  equivalent,  as  shown  elsewhere.^ 


V.  ABSORPTION  PROFILE 

From  Eq.  (1.1)  we  notice  that  the  absorption  profile  involves  quantities  as 
Trp(t)/i(t+r )n(t) ,  which  depend  on  two  times,  and  their  evaluation  is 
correspondingly  more  complicated.  In  this  section  we  derive  a  formal  expression 
for  I(w). 

First  we  introduce  a  two-time  operator 


D(r,t)  -  e'lLr(MlP(t) 


(5.1) 


and  its  Fourier-Laplace  transform 


D(w, t) 


dr  elwr  D(r , t)  , 


(5.2) 


where  the  frequency  dependence  only  refers  to  the  first  argument  of  D(r,t).  We 
remark  that  D(r,t)  is  an  operator  in  the  entire  atom  plus  crystal  Hilbert  space. 
Comparison  of  Eqs .  (5.1)  and  (5.2)  with  Eq.  (1.2)  shows  that  the  expression  for 
I(w)  can  be  cast  in  the  form 


I(u)  -  I_ | e  • « 

L.  —  Z  — 


L'  « 


Re  lim  Tr^D(w, t) 
t"*« 


(5.3) 


From  the  fact  that  the  dipole  operator  p  acts  only  on  adsorbate  states,  it  follows 
chat 


TrpD(w.t)  -  Tr  oD  (w.t)  ,  (5. A) 

a  a 

wi-ch  Da(w,c)  —  TrpD(w,t),  the  reduced  adsorbate  operator. 

Then  we  notice  that  D(r,t)  obeys 

i  ^  D(r,t)  -  L  D(r,t)  ,  r  >  0  ,  (5.5) 

as  is  evident  from  its  definition  (5.1).  Hence  the  r -dependence  of  D(r,t)  is 
governed  by  the  same  equation  as  the  one  which  determines  the  time  evolution  of 
the  density  operator  p(t)  ,  Eq.  (1.4).  An  important  difference  is  that  the  initial 
value  for  Eq.  (5.5)  reads 

D(0 , t)  -  [p,p( t)]  ,  (5.6) 

in  contrast  to  the  equation  for  p( t) ,  where  p( 0)  can  be  chosen  arbitrarily.  The 

time  regression  of  D(r,t)  on  the  interval  0  <  r  <  ®  is  identical  to  the  time 

evolution  of  p(t)  on  0  <  t  <  ®,  and  therefore  it  should  be  possible  to  express 

D(w,t)  in  r(w) .  Additionally,  the  initial  value  D(0,t)  depends  explicitly  on  t, 

via  p(t),  and  in  a  finite  memory-time  theory  this  quantity  will  carry  a  memory  to 

the  time  evolution  of  p(t)  in  the  recent  past  (times  smaller  than  t)  . 

Contributions  to  the  line  profile  which  arise  due  to  this  mechanism  will  be 

referred  to  as  initial  correlations.  It  might  seem  that  in  the  limit  t  -*  ®  we  can 

replace  p(t)  by  p  p  ,  which  would  eliminate  initial  correlations.  We  shall  show 
a  p 

that  this  is  not  correct  in  general.  In  Ref.  15  we  have  developed  a  general 


1 

method  for  the  evaluation  of  quantities  of  the  form  D  If  we  apply  that 

3. 

theory  to  the  present  situation,  we  find  the  formal  expression  for  the  absorption 
profile 


I(B)  “  II'  — z*— Ll 


"jte  ReTra4  <ViT<«»'« 


(5.7) 


where  the  Liouvillians  L,,  and 


and  Ly 


Va  “  a  ’ 
Va  "  [^’aa] 


(5.8) 

(5.9) 


are  introduced  in  order  to  simplify  the  notation.  Equation  (5.7)  is  the  most 
condensed  and  general  expression  for  the  absorption  profile  of  an  atomic  bond  on  a 
crystal.  We  recognize  the  time  regression  operator  i/(w-L  +ir(w)),  which  is  the 

3 

same  indeed  as  in  Eq.  (4.1).  Initial  correlations  are  accounted  for  by  the 
frequency-dependent  operator  T(w),  which  is  explicitly 


T(u)c7  -  Tr  L  ~  j 
a  p  ap  w-L  -L 


t  —1 -  L  (op) 

^  iO+-L  -L  aP  a  P 


(5.10) 


Furthermore,  we  notice  that  the  series  of  Liouvillians  under  the  trace  in  Eq. 
(5.7)  act  on  the  stead-state  density  operator  p a  of  the  adsorbate,  which  can  be 
obtained  by  solving  Eq.  (4- A). 


VI.  MATRIX  ELEMENTS  OF  T(d)  AND 


Eigenvalues  and  eigenstates  of  the  atomic-bond  Hamiltonian  H  from  Eq.  (2.2) 

cl 

are  defined  by 

Hjp>  -  tf«p|p>  .  (6.1) 

and  due  to  the  neglect  of  lateral  motion  the  eigenvalues  are  non- degenerate . 
For  realistic  adsorbate  systems  there  are  approximately  25  bound  states  |p>.  In 
this  section  we  expand  the  various  Liouvillians  onto  the  set  { | p> )  of  adsorbate 
bound  states. 

Without  coupling  to  the  reservoir  (T(w)  -  0,  T(u>)  -  0)  the  time  regression 
operator  is  determined  by  the  inverse  of  w-L  ,  and  with  Eq .  (6.1)  we  readily  find 

<p|  {(w-LJaj  |q>  -  (cj-A  )  <p  I cr I q>  ,  (6.2) 

a  pq 

in  terns  of  the  level  separations 


(6.3) 


Equation  (6.2)  relates  the  matrix  elements  of  (u-L^)a  to  the  matrix  elements  of  a 
for  any  adsorbate  density  operator  a,  and  thus  Eq.  (6.2)  implies  the  matrix 
representation  of  the  Liouvillian  w-La< 

The  coupling  between  the  crystal-atom  motion  and  the  adsorbate  motion  is 
established  by  the  Hamiltonian  H  of  Eq.  (2.3),  which  has  the  adbond  part 


A  ‘  /  V  lA.I 


Matrix  elements  of  this  Hermitian  operator  will  be  denoted  by 


S  -  <p  S  q>  -  S* 
pq  vi  i  qp 


(6.5) 


Similarly  we  denote  the  matrix  elements  of  /i  by  ^pq  -  <p|/i|q>  —  /i*p .  case  °f  a 

Morse  potential  explicit  expressions  for  Spq  and  /jpq  can  be  derived,  with  an 
integral  due  to  Rosen. ^  Furthermore,  we  remark  that  for  any  potential  V(z)  the 
diagonal  matrix  elements  of  its  derivative  vanish  with  respect  to  the  eigenstates 
of  Ha,  e.g. ,17 


S  -  0 
PP 


(6.6) 


The  permanent  dipole  moments  pipp  of  level  |p>,  however,  are  finite  in  general. 

With  the  methods  of  Ref.  15  we  can  evaluate  the  matrix  elements  of  T(w) ,  as 
defined  in  Eq.  (4.2),  with  respect  to  the  basis  set  (|p>).  The  results  is 

<p| (r(w)tr) |q>  -  ^  (?(Aqa+«)SpaSab<b|a|q>  +  l*(Apa-w)S*aS*b<p|a|b> 


‘  ^V^Vbq  +  ?*<Aaq-w>SVaP)<a|ff|b>1 


(6.7) 


in  terms  of  the  matrix  elements  of  S  and  the  reservoir  correlation  function  5(w) 
from  Section  III.  It  appears  that  the  frequency  dependence  of  T(w)  enters  as  a 
shift  of  the  level  separations  Apq  in  the  arguments  of  the  correlation  function. 

The  initial  correlation  operator  T(w)  from  Eq.  (5.10)  involves  the  inversion 
of  two  Liouvillians .  After  laborous  computations  we  obtain 


<P |  {T(<j)<7 )  | q>  -  )  )-?(Aqb+o,))SpbSca<a|a|q> 


-  (?*(ACQ)  -  ?*(\q-))SpbSaq<ci<7ia> 


+  (?(Abp)-?(Acp+W))ScqSpa<a|a|b> 


'  (?*(Apb)-^(V-w))ScqSab<plfl!a>) 


— 2^-  l(f(i  )  -  f  (A  +a>))S  S  ,  <b  |  cr  a> 
A  +w  v  ac'  ap  aq  cb  11 


-  <I*<Aca)  ‘  J*(Apa'")  ’  SaqSba<C  I "  I b>  1 


\b*“ 


tf<\a>  •  I(Aq.4“))Sp.S«c<ck|b> 


(i*(Aab)  -  J*(4aq..a))SpaScb<a|a|c>} 


(6.8) 


Again,  the  frequency  depenence  enters  as  a  shift  of  the  arguments  in  ?(Aa^) ,  but 

in  addition  overall  factors  (A  ,  +w)  ^  appear.  Furthermore,  we  see  that  T(w) 
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depends  on  the  matrix  elements  of  the  dipole  operator. 


With  the  matrix  representations  of  the  various  Liouvillians  it  is 

straightforward  to  evaluate  I(w)  from  Eq.  (5.7)  for  any  configuration  of  levels, 

or,  for  any  potential  V(z).  The  profile  I(w)  will  exhibit  many  overlapping  lines 

at  the  adsorbate  resonances  A  ,  .  In  order  to  disentangle  the  contributions  to 

ab 


I(w)  from  the  different  transitions,  and  to  elucidate  the  significance  of  initial 
correlations,  we  elaborate  on  the  situation  where  the  potential  supports  only  two 
bound  states.  Let  us  denote  these  states  by  |2>,  |l>,  with  the  convention  that 


1 


'->0  w2  *  W1  >  0-  from  Eq.  (6.5)  it  follows  that  -  S ^  ”  0,  and  therefore  the 


only  non-vanishing  matrix  element  of  dV/dz  is 


SQ  “  <1 1 S | 2>  , 


(7.1) 


which  will  be  assumed  to  be  real  (as  for  a  Morse  potential). 

From  Eqs .  (6.2)  and  (6.7)  we  find  the  matrix  representation  of  w-L  +iT(w) . 

cl 

On  the  basis  1 2x2 1  ,  |lxlj,  |2xl|,  1 1x2 1  this  becomes 


w-L  +  ir(w)  - 

cL 


ort-ip(w)  - iq(w) 

-ip(w)  aH-iq(w) 


w-w  +in(w) 
o 

-  ir?(w) 


-ir7(w) 


w+a>o+iry  (w) 


(7.2) 


in  terms  of  the  parameter  functions 


p(w)  -  Sq  (?(«o+«)  +  ?*(wq-w)) 


(7.3) 


q(w)  -  Sq  (?( -wq+w)  +  ?*(-wq-w)) 


(7.4) 


»j(w)  -  SqZ(?(w)  +  ?*(-w)) 


(7.5) 


Inversion  of  the  matrix  (7.2)  then  yields  the  resolvent  for  fi a(u) .  Eq ■  (^D.  anci 
a  Fourier -Laplace  inverse  of  the  result  gives  pa(t) •  Recalling  the  complicated 
frequency  dependence  of  ?(w),  Section  III,  shows  that  an  evaluation  of  pfl(t)  in 
the  transient  regime  0  <  t  <  ®  is  evidently  impossible,  unless  numerical  methods 
are  applied. 

Fortunately,  the  absorption  profile  depends  only  on  the  steady-state  density 
operator  p  (t  —  <*>) ,  which  obeys  Eq.  (4.4).  With  the  matrix  representation  (7.2) 
this  equation  is  easily  solved,  with  result 


p a  -  | 2>n^<2 |  +  |l>n^<l| 


(7.6) 


The  steady-state  level  populations  are 

_  - alO) -  _  - ZlQl -  (7.7) 

°2  P(0)  +  q(0)  '  nl  p«3)  +  q(0)  ’ 

and  Eq.  (7.6)  expresses  that  the  coherence  <l|p  | 2>  vanishes,  as  usual  in  thermal 

3. 

equilibrium.  From  the  definitions  of  p(0)  and  q(0)  in  Eqs.  (7.3)  and  (7.4)  we  see 
that  the  factor  Sq  drops  out,  and  hence  the  populations  are  completely  determined 
by  the  reservoir  correlation  function  ?(w),  at  the  resonance  w  -  wQ-  With  Eq. 
(3.6)  we  then  find  that  the  dependence  on  £  also  disappears,  so  that  nj  and  are 
determined  by  g(z) ,  which  has  only  the  temperature  as  parameter.  Explicitly, 


_ Re  g(-(dQA;D) _ 

n2  "  Re(g(wQ/wD)  +  g(-wQ/wD)) 

Re  g(wQ/wn) _ 

nl  “  Re(g(u>  /w  )  +  g(-w  /w  )) 


(7.8) 


(7.9) 
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which  contain  only  the  real  parts  of  the  correlation  functions.  Then  it  follows 


from  Eq .  (3.10)  that  for  the  level  populations  are  undetermined.  This  is  a 


consequence  of  the  fact  that  we  restricted  the  atom-crystal  interaction  to  single¬ 


phonon  couplings. 


VIII. 


For  a  two-state  system  the  initial  correlation  operator  T(w)  is  a  4  x  4 


matrix,  but  its  general  representation,  as  it  follows  from  Eq.  (6.8),  is  still 


cumbersome.  For  the  evaluation  of  the  spectrum  I(u),  however,  we  only  need  to 


know  the  result  of  its  action  on  p  ,  as  can  be  seen  from  Eq .  (5.7).  Since  p  has 
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only  two  non- vanishing  matrix  elements,  rather  than  four,  this  simplifies  the 


situation.  On  the  same  basis  as  the  representation  (7.2),  we  then  find 


T(w)pa  - 


(8.1) 


-y(w) 


x(w)  -  —  (mu-m22) 


x  {n»(?(w  )  -  ?(w  Hi>)  -  l*(u>  )  +  ?*(u>  -u>)) 
Zoo  o  o 


+  n.(?(-u>  )  -  ?(-u>  +w)  -  ?+( -w  )  +  ?*(-CJ  -U)))) 
loo  o  o 


(8.2) 


y(w)  - 


■m  »-r  r*  vr  -mj*  tt-»  -mj~,  w^-\  ^ 


2  2  m21 
O)  -CJ 

o 


x  (n2((?(wQ)  -  ?(w) ) (wq-kj)  +  (f*(wQ)  -  I*( -u' ) (w-wq) ) 


+  n^((l(-wQ)  -  ?((j))(wq-w)  +  (?*(-w)  -  f*( -wq)  )  (w-hjq)  )  ) 


(8.3) 


The  right-most  factor  in  Eq.  (5.7)  is  (Ly- iT(w) )p^  ,  and  therefore  the 
relative  significance  of  the  initial  correlations  follows  from  a  comparison  of 
T(w)^a  with  Lyp^.  With  Eq.  (5.9)  we  find 


Va  “ 


(8.4) 


The  most  important  difference  is  that  IyPa  only  depends  on  the  transition  dipole 

matrix  element  whereas  T(w)  acquires  a  contribution  from  the  permanent  dipole 

moments  of  the  two  levels  (terms  proportional  to  x(w)).  Furthermore,  T(o>)  is 
2 

proportional  to  Sq  (the  strength  of  the  interaction  between  the  vibrating  atom 
and  the  crystal),  whereas  IyPfl  is  independent  of  this  parameter. 


An  absorption  profiel  for  a  two- state  system  is  called  a  line  shape,  because 
it  singles  out  a  specific  transition  of  the  vibrational  bond.  With  the  matrix 
representations  of  the  previous  sections  we  are  now  able  to  construct  the  line 
shape  I(u)  by  simple  matrix  operations.  First  we  define  the  dimensionless  line 
shape  I(u)  by 


-j.  n 


I(w)  -  I(w){I  (fo)<c)  l£z,lLf  ^21' 


in  order  to  suppress  irrelevant  overall  factors.  Then  we  write  I(w)  as  a  sum  of 
two  contributions 


x<«>  -  1^“>reg  +  I(w)in  * 


(9.2) 


where  I(w)^e^  comes  from  the  term  (regression  part)  and  I(w)^n  represent  the 
initial-correlation  contribution.  It  appears  that  I(w)  can  be  expressed  entirely 
in  the  parameters 


a  -  M 


a  „  <*22^2 

"2! 


(9.3) 


(9.4) 


Here,  a  equals  the  half  width  at  half  maximum  if  the  line  would  be  approximated  by 
a  Lorentzian,  and  in  the  case  T  -  0,  wq  -  and  m  measures  the  relative 

importance  of  the  permanent  dipole  moments  in  comparison  with  the  transition 
dipole  moment.  We  finally  obtain 


I(w)  -  2 uu  (n,  -n„)  Re  ~ — f - 

'reg  ov  1  2'  D(w>coh 


I(w) in  -  2wwo  Re,D(u>)  7  *(w)  '  D(w)  *(w)) 

coh  pop 


(9.5) 


(9.6) 


in  terms  of  the  auxiliary  functions 


D(u)coh  “  w  '  Wo  +  - 


(9.7) 


D(«)pop  ”  w(oH-i(p(w)+q(w)))  , 


(9.8) 


g("n)-g(“)  g*(-w)-g*(w  ) 

-  -2ia{iu( -  +  - — - —) 

2.  u  ~u  w+w 

o  o 


g(-a>o)-g(w)  g*(-UQ)-g*(-w) 

+  ni^  Ci)  -Hi)  U)  -Ci)  )  )  » 

o  o 


(9.9) 


tf(«)  “  2^”  (n2(g(“o)  '  S(«0+«)  ‘  S*^0}  +  g*("0-^)) 


+  n1(g(-<jQ)  -  g(-wo+o>)  -  g*(-wQ)  +  g*(-£>o-<i)) } 


(9.10) 


with  C>  -  W/WD  the  frequency  in  units  of  ci^.  The  two  D(cj)  functions  are  the 
determinants  of  the  two  submatrices  in  Eq.  (7.2),  and  the  subscripts  coh  and  pop 
refer  to  the  subspaces  of  coherences  ( |  2x1 1  ,  1 1X2  | )  and  populations  (  |  2x2  |  , 
|lxl|),  respectively.  The  above  sequence  of  formulas  determines  the  shape  of  a 
single  absorption  line,  where  wq,  a,  7  and  m  are  the  only  arbitrary  parameters  (if 
we  take  as  the  frequency  unit) .  For  a  specified  potential  V(z)  we  can  express 
cjQ,a  and  ra  in  properties  of  this  potential,  and  7  is  simply  the  temperature  of  the 
crystal. 

Substitution  of  ^(w)co^  into  Eq.  (9.5)  and  removal  of  the  real  part  yields 


4 ( n^  -  ) cjqu>  RerKu>) 


rC®  (w  ^ -ci)^+2(i)Imr; (w)  )  ^  +  4w^  (Ren (w) )  ‘ 


(9.11) 


From  the  definition  (7.5)  of  r/( w)  and  the  property  (3.10)  of  the  reservoir 
correlation  function,  it  follows  that  Re>?(t*>)  -  0  for  w  >  and  consequently 


I(w>reg  “  0  ’  for  w  >  WD 


(9.12) 


Therefore,  absorption  for  u  >  can  only  be  a  result  of  nonvanishing  initial 
correlations . 


Before  we  discuss  the  relevance  of  the  finite  memory- time  in  the  time 
regression  of  dipole  correlation  functions  and  the  importance  of  the  inclusion  of 
initial  correlations,  we  summarize  the  results  from  earlier  theories.  A  most 


obvious  approximation  would  be  the  factorization 


lim  D(O.t)  „  [/i,p  l/>  (10.1) 

t-«o  3  p 

of  the  initial  condition  for  the  time-regression  equation  (5.5).  It  can  be 
shown^  that  this  implies  T(w)  —  0,  and  consequently  I(w).^  ”  0,  so  that  I(w)  is 
approximated  by  I(w)reg. 

A  more  rigorous  simplification  is  the  Markov  approximation,  in  which  any 
memory  effect  is  discarded.  First  this  implies  the  factorization  (10.1),  where  a 
memory  of  p( t)  to  its  recent  past  is  neglected.  Secondly,  we  adopt  a  memoryless 
description  of  the  time  regression,  which  yields  a  frequency- independent 
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relaxation  operator  f  .  The  Markovian  equivalent  of  the  matrix  (7.2)  reads 


a 


m 


uH-ip(O) 
- ip(O) 


• iq(O) 
ott-iq(O) 


cj-o>0+if)  (WQ) 

-i» ?<w0) 


-ifj*(wo) 

urt-w  +  ir7*(<w  ) 
o  o' 


(10.2) 


In  the  left-top  submatrix,  the  parameter  functions  p(w)  and  q(w)  are  replaced  by 
their  values  at  w  *  0,  Since  the  steady-state  level  populations  n^  and  nj  are 
determined  by  T(0) ,  we  conclude  that  n^  and  ^  are  unaffected  by  this 
approximation.  In  the  right-bottom  submatrix  the  functions  rj(u)  are  now  evaluated 
at  the  resonance  frequency  ,  rather  than  at  the  laser  frequency  oj,  and 
furthermore,  the  functions  r/(wQ)  in  the  fourth  column  are  now  complex-conjugated. 
In  this  approximation  the  line  shape  is  found  to  be 


I(W)  - 


4(n^-n2)wQw  Rer7(u>Q) 


m  (w  ^-td^+2w  Imr;(w  ))^  +  4w  (Ref?(w  )) 


(10.3) 


which  greatly  resembles  I(u>)  from  Eq.  (9.11).  The  most  important  difference  is 

reg 

that  we  now  find  Rerj(wQ)  in  the  numerator,  rather  than  Rery(w)  ,  which  implies  that 
I(w)  is  finite  for  w  >  wn. 

in  u  i 

In  the  most  simple  theory  of  relaxation  (which  leads  to  a  master  equation  for 

the  time  evolution  of  the  populations),  we  furthermore  neglect  any  coupling 

between  coherences  and  populations,  and  between  coherences  which  evolve  with  a 

different  frequency  in  a  free  evolution  (no  coupling  to  the  reservoir) .  This 

approximation  is  usually  called  the  secular  approximation,  and  it  is  equivalent 

to  the  neglect  of  terms  which  oscillate  fast  on  a  time  scale  1/a  (Eq.  (9.3)). 


2- 


I 


•-w  "Ory  »vp  ^Tipr^T^rv*:  w '  nr.-nr-'  »»>r  m-M  wt  w  w~  »v  wv  -w--. 


Then  we  replace  the  off-diagonal  matrix  elements  -i  »?(wo)  anc*  -i>7*(w0)  in  Eq . 
(10.2)  by  zero,  which  gives  for  the  line  shape  (s  -  secular) 


I(«), 


w  (n.  -n„)  Re  - — — — r 

o  1  2  d-(j)  +ifj(w  ) 


(10.4) 


We  obtain  the  standard  Lorentzian,  which  has  its  maximum  at  w  -  coq  +  Imi7(w  ),  and 
has  a  half-width  at  half-maximum  equal  to  RerKwQ) • 

It  is  easy  to  verify  that  the  three  approximations  yield  the  same  value  for 
I(u>)  at  the  resonance  frequency  ,  e.g., 


I(w  )  -  I(w  )  -  I(w  ) 

o  reg  o  m  os 


wo(ni-n2)Re  ^TT 


(10.5) 


Conversely,  this  implies  that  the  successive  refinements  will  have  a  major 
significance  in  the  line  wings  only,  which  is  illustrated  in  Fig.  1. 

XI.  RESULTS 

The  line  shape  I(w)  is  composed  of  two  contributions,  which  are  drawn 

separately  in  Fig.  2  for  a  specific  set  of  parameters.  We  notice  that  the  initial 

correlation  part  is  not  small,  and  that  it  takes  on  both  positive  and  negative 

values.  A  considerable  increase  of  the  absorption  is  found  in  the  line  center  and 

the  red  wing,  whereas  the  blue  wing  is  only  slightly  modified,  both  in  comparison 

with  the  line  shape  I(w)  (the  best  approximation  so  far).  Due  to  the  neglect 

reg 

of  permanent  dipole  moments  in  this  case  (m  -  0) ,  we  find  a  sharp  edge  at  u>  -  u>^, 
and  a  vanishing  absorption  for  w  >  oi^.  Of  course,  multiphonon  processes  would 
also  give  rise  to  absorption  at  u  >  but  these  contributions  are  assumed  to  be 
small.  The  cutoff  can  be  understood  from  the  fact  that  a  photon  absorption  from 
the  laser  gives  rise  to  a  transition  | 1>  -*  ) 2>  of  the  adsorbate.  Energy 
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conservation  then  implies  that  this  process  must  be  followed  by  a  single -phonon 
emission  into  the  crystal,  under  a  decay  j 2>  -»  j 1>  of  the  adsorbate  bond.  In  a 
Debye  model  there  are  no  phonon  modes  available  for  w  >  and  therefore  this 
process  cannot  occur. 

A  remarkable  profile  arises  if  we  allow  the  atomic  bond  to  have  a  permanent 

dipole  moment  (m  *  0) ,  which  is  illustrated  in  Fig.  3.  For  u  <  the  value  of  m 

has  hardly  any  significance,  but  above  the  cutoff  frequency  we  now  find  a  finite 

absorption  if  m  **  0.  We  observe  a  peak  at  co  —  a  smooth  background  for 

co  <  u>  <  u  4-  t .  How  can  this  be  understood?  First  we  remember  that  a  low- 
D  o  D 

intensity  profile  is  a  balance  between  the  stimulated  absorption  and  emission 
rates  for  single-photon  transitions.  Hence  the  phenomenon  cannot  be  attributed  to 
multiphoton  processes.  Then  we  recall  that  we  restricted  ourselves  to  a  model  of 
single-phonon  interactions,  which  rules  out  multiphonon  processes.  Third,  a 
positive  I(tj)  corresponds  to  a  photon  absorption.  A  positive  I(w)  in  the  range 
w  >  >  u  can  then  only  be  found  from  the  following  energy-conserving  process. 

Initially  the  adsorbate  is  in  its  lower  state  |l>-  Absorption  of  a  photon  then 
excites  the  system  to  a  virtual  level  with  energy  J’fw  above  the  energy  of  state 
|l>,  and  subsequently  this  state  decays  to  | 2> ,  accompanied  by  the  emission  of  a 
phonon  with  frequency  w  -  w  .  The  second  transition  can  only  occur  if  u>  -  u>q  < 

Up,  or  «  <  u>q  +  Wp.  This  explains  why  I(w)  vanishes  identically  for  o>  >  u>q  + 

I 

and  can  be  finite  for  con  <  w  <  w  +  wn. 

D  o  u 


XII.  CONCLUSIONS 

We  performed  very  detailed  calculations  on  the  optical  spectral  profile  and 
line  shapes  of  physisorbed  atoms  on  the  surface  of  a  harmonic  crystal.  Coupling 
of  the  atomic  motion  to  the  phonon  field  of  the  crystal  provides  the  relaxation 
mechanism  for  the  evolution  of  the  adsorbate  vibrational  states  towards  thermal 
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equilibrium.  Since  the  amplitude  correlation  function  of  a  substrate  atom  has  a 
finite  decay  time,  the  time  evolution  of  the  reduced  adsorbate  density  operator 
exhibits  a  memory  effect.  The  absorption  profile  is  determined  by  dipole 
correlation  functions,  depending  on  two  times.  We  identified  two  different 
aspects  of  the  memory,  which  were  called  regression  and  initial  correlation.  A 
general  finite  memory- time  reservoir  theory  was  applied  for  the  evaluation  of  the 
line  shape. 

It  appeared  that  the  properties  of  the  reservoir  could  be  accounted  for  by  a 
single  dimensionless  function  g(z) ,  which  has  only  the  dimensionless  temperature  7 
as  parameter  (apart  from  a  scaling  factor  £) .  The  other  parameters  of  the  line 
shape  are  the  resonance  frequency  of  the  uncoupled  adsorbate,  the  interaction- 
strength  parameter  a  (frequency,  related  to  line  width),  and  the  permanent-dipole 
parameter  m  (dimensionless).  From  our  analytical  expression  for  the  line  shape  we 
showed  that  any  absorption  above  the  Debye  frequency  can  only  be  a  consequence  of 
nonvanishing  initial  correlations  (in  the  single-phonon  approximation).  It  was 
exemplified  (Fig.  3)  that  I(w)  is  finite  indeed  for  u  >  provided  that  the 

transition  has  a  permanent  dipole  moment  (as  is  the  case  for  a  Morse  potential) . 
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Fig .  1 .  Secular  (a),  Markov  (b)  and  regression  (c)  approximacion  of  I(w)  as  a 
function  of  £  -  and  for  -  0.65  x  (dotted  line),  a  -  0.4  x  7  -  10 

(low  temperature)  and  m  -  0.  Curve  a  is  a  Lorentzian  around  +  ImrKw^)  ,  but  for 
the  present  parameters  the  line  shift  is  negligible  (Imrf(w  )  -  -0.011).  As 
predicted  in  Eq.  (10.5),  the  value  of  I(u)q)  is  the  same  for  the  three  curves.  The 
secular  approximation  gives  a  finite  absorption  for  u  0,  whereas  the  improved 


results,  curves  b  and  c,  give  a  vanishing  absorption  for  <j  -*  0 .  In  the  blue  wing, 
however,  curves  a  and  b  remain  finite,  but  the  most  refined  theorv  gives  I(~>)  -  0 
for  w  >  cjp. 

Fig .  2 .  Curve  a  and  b  are  I(w)  and  I(w).  ,  respectively,  for  the  same 

reg  in  j 

parameters  as  in  Fig.  1.  Their  sum,  curve  c,  is  I(u>),  which  is  calculated  with 
the  present  theory.  Comparison  of  curves  a  and  c  shows  the  significance  of  the 
improvement.  Note  also  the  considerable  red  shift  of  the  line  with  respect  to  thi 
resonance  frequency  . 

Fig,  3 .  Curve  a  is  the  same  as  curve  c  from  Fig.  2  (m  -  0) ,  and  for  curve  b  we 
took  the  dipole -moment  parameter  ro  equal  to  0.7.  The  nonzero  value  of  m  appears 
to  have  no  importance  at  all  for  u  <  but  for  w  >  only  the  profile  b  remain! 

finite.  As  pointed  out  in  the  text,  the  absorption  for  w  >  is  entirely  due  to 


the  initial  correlations. 
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